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Abstract: We study a model of "organized" criticality, where a single avalanche propagates 
through an a priori static (i.e., organized) sandpile configuration. The latter is chosen accord- 
ing to an i.i.d. distribution from a Borel probability measure p on [0, 1]. The avalanche dynamics 
is driven by a standard toppling rule, however, we simplify the geometry by placing the problem 
on a directed, rooted tree. As our main result, we characterize which p are critical in the sense that 
they do not admit an infinite avalanche but exhibit a power-law decay of avalanche sizes. Our anal- 
ysis reveals close connections to directed site-percolation, both in the characterization of criticality 
and in the values of the critical exponents. 
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1. Introduction 

1.1 Motivation. 

Since its discovery by Bak, Tang and Wisenfeld [1, 2], self-organized criticality (SOC) has re- 
ceived massive attention in the physics literature. Variants of the original sandpile model of [1] 
were studied and some of them even "exactly" solved (see [4, 5, 10, 1 1] or [7] for a recent review 
of the subject). However, despite great efforts and literally thousands of published papers, the 
present mathematical understanding of SOC lags far behind the bold claims made by physicists. 
Much of that failure can be attributed to the fact that the models used to demonstrate SOC are 
difficult to formulate precisely and/or too difficult to study using the current techniques of prob- 
ability theory and mathematical physics. From the perspective of the latter fields, the situation 
seems ripe for considering models which concern at least some aspects of SOC, provided there is 
a decent prospect of a self-contained rigorous analysis. 

The general idea behind SOC models is very appealing. Consider for instance Zhang's sandpile 
model [12] on Z 2 , where each site has an energy variable which evolves in discrete time-steps 
according to a simple "toppling" rule: If a variable exceeds a threshold value, the excess is 
distributed equally among the neighbors. The neighboring sites may thus turn supercritical and 
the process continues until the excess is "thrown overboard" at the system boundary. What makes 
this dynamical rule intriguing is that, if the toppling is initiated from a "highly excited" state, 
then the terminal state (i.e., the state where the toppling stops) is not the most stable state, but 
one of many least-stable, stable states. Moreover, the latter state is critical in the sense that 
further insertion of a small excess typically leads to further large-scale events. Using the sandpile 
analogy, such events are referred to as avalanches. 

In this paper, we study the scaling properties of a single avalanche caused by an overflow at 
some site of a critical (i.e., least-stable) state. However, as indicated above, the full problem is 
way too hard and we have to resort to simplifications. Our simplifications are twofold: First, we 
treat the energy variables of the critical state as independent and, second, we consider the model 
on a directed, rooted tree rather than 1? . The first assumption is fairly reasonable, at least on a 
coarse-grained scale, because numerical results [6] suggest a rather fast decay of spatial correla- 
tions in the critical states. The second assumption will allow us to treat the correlations between 
different branches of the avalanche as conditionally independent, which will greatly facilitate the 
analysis. Finally, the reduced geometry allows for the existence of a natural monotonicity not 
apparent in the full-fledged model. 

While placing the model on a tree simplifies the underlying geometry, some complexity is 
retained due to the generality of the single-site energy variable distribution. In fact, the set of 
underlying distributions plays the role of a parameter space in our case. As our main result, 
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we characterize the subspace of distributions for which the configurations of energy variables 
have exactly the behavior expected from the SOC states: no infinite avalanches but a power- 
law decay of avalanche sizes. As it turns out, there is a close connection to site-percolation on 
the underlying graph, both in the characterization of criticality and in the values of the critical 
exponents. However, the significance of this connection for the general SOC models has not yet 
been clarified. 

1.2 The model. 

In order to precisely define our single-avalanche model, we need to introduce some notation. 
Let b > 1 be an integer and let T5 be a &-nary rooted tree, with the root vertex denoted by 0. 
We use \a\ = k to denote that a G T& is on the k-th layer. When \a\ = k, we represent a as 
a fc-component object. Each component is an integer in {1, . . . , b}; hence the site label can be 
used to trace the path from a back to the root. If a is an £-th level site with I > 0, we let m(a) 
denote the "mother-site." Explicitly, if a = (a±, . . . , og), then m(a) = (01, ... , cr^-i). The 
edges of T& are the usual directed edges { (V, a) G T& x T& : a' = m(cr)}. 

Let M. be the space of all probability measures on the Borel u-algebra of [0, 1]. Fix a p G M. 
and let F p = p Tb . Let E p denote the expectation with respect to P p . The dynamical rule driving 
the evolution is defined as follows: Let X = (X a ) ae j b be the collection of i.i.d. random variables 
with joint probability distribution P p and let v G (0, 00). The process generates the sequence 

X^(t) = (X^(t)) aeJb , 4 = 0,1,..., (1.1) 
obtained from the initial condition 

X^(O) = l X0 + V ' ifa = ' (1-2) 

I X a , otherwise, 

by successive applications of the deterministic (Markov) update rule 



X^{t + l) = { 



[X£\t) + \X^ (t), ifX m{a) (t)>l, 



0, ifxj?\t)>l, (1-3) 

K Xa\t), otherwise. 



Note that, if X^ v \t + 1) = X^\t) for all a G T b , then X^\t) < 1 and the process has 

effectively stopped. (However, we let xj?\i) be defined by (1.3) for all t > 0.) 

Here is an informal description of the above process: Starting at the root we first check 
whether X + v > 1 or not. If not, the process stops but if so, then this value is distributed evenly 
among the "daughter" cells, which have their values updated to xj?\l) = X a + \{X + v). 

The value X^\l) is set to zero and we say that the root has "avalanched." If none of the up- 
dated "daughter" values exceed one, the process terminates; however, if there is any first-level a 
with xj?\l) > 1, then X^\l) is set to zero, the value xj?\l) is evenly distributed among the 
"daughters" of a and we say that a has "avalanched." The process at future times is described 
similarly. 
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Obviously, the variables X a play the role of the "energy variables" in the description of 
Zhang's avalanche model in Section 1.1. In our case the critical threshold is one, but, in (1.3), 
we chose to distribute the entire value of an "avalanching" site rather than just the excess to the 
(forward) neighbors. This choice is slightly more advantageous technically. 

1.3 Main questions and outline. 

Let s/W{t) = {a € T b : \t) = 0, X^ v \s) / for some s < t} be the set of sites that 
have "avalanched" by time t. Similarly, let srf^ = IJt>o W ^ e tne set °f s i tes tnat wu ^ ever 
avalanche. We use l^^l to denote the number of sites in the avalanched set (which includes the 
possibility of \srf^ | = oo). The set srfW and its dependence on p and v are the primary focus of 
our study. 

The first question is whether the process X^ (t) lives forever, i.e., is there an infinite avalanche? 
More precisely, for what measures p G M. is the probability 

A&>=P p (|^M| = oo) (1.4) 

non-zero for some value of vl A related question is whether the average size of the avalanched 
set is finite. The relevant object is defined by 

X ^=E p (|^)|). (1.5) 

(Notice that, due to the directed nature of the dynamical rule, both quantities and x^ w& 
monotone in the underlying measure and v.) Again, we ask: For what measures p we have = 
oo for some v ? In addition, we might ask: Is the divergence of the mean avalanche size equivalent 
to the onset of infinite avalanches or can there be an intermediate phase? 

To give answers to the above questions, we will parametrize the set M by values of a partic- 
ular functional 3: M. — ► [0, 1]. Here i(p) roughly corresponds to the conditional probability in 
distribution P p that, given the avalanche has reached a site a G T5 far away from the root, the 
site a will also avalanche. (The definition of 3 is somewhat technical and we refer the reader 
to Section 2.2 for more details.) The characterization of the avalanche regime in terms of 3 is 
then very transparent: There is a critical value } c = \, such that the quantity x^ f° r measure p 
diverges if i(p) > 3 C and v is sufficiently large, while it is finite for all v if i(p) < 3 C . Sim- 
ilarly we show, for a reduced class of measures, that for measure p vanishes for all v if 
and only if i(p) < % c . These results are formulated as Theorems 2.4 and 3.1 in Sections 2.2 
and 3.1, respectively. (Outside the reduced class of measures, there are some exceptions to the 
rule that = for measures p with }{p) = 3 C , i.e., there are some measures which avalanche 
also at criticality, see Remarks 1 and 2 in Section 2 for more details. These examples are fairly 
contrived, so we exclude them from further considerations.) 

Note that, for both quantities (1.4) and (1.5), the transitions happen at the same value, 3 C , 
which rules out the possibility of an intermediate phase. To elucidate the behavior of 3 near 3 C , it 
is worthwhile to introduce appropriate critical exponents. In particular, we ask whether there is a 
critical exponent 7 > such that 

X [v) ~(3c-3(p)r 7 , 3(p)t3c (1-6) 
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an exponent [3 > such that 

^ ) ~(3(/>)-3c) /3 , 3(p)l3c (1-7) 
and, finally, an exponent 5 > such that if %{p) = j c , then 

¥ p {W {v) \ > n) ~n- 1/5 , n^oo. (1.8) 

All of these three relations of course include an appropriate interpretation of the symbol "~" and, 
with the exception of the last relation, also an interpretation of the limit "3(p) tends to 3 C ." 

The relations for the critical exponents are the subject of Theorem 4.1 in Section 4. The upshot 
is that all three exponents take the mean-field percolation values, 

7 = 1, (3 = 1, 5 = 2. (1.9) 

Neither the fact that the critical value 3 C equals the percolation threshold for site percolation on T& 
is a coincidence. Indeed, the avalanche problem can be characterized in terms of a correlated- 
percolation problem on T& (see Section 2). 

We finish with a brief outline of the paper: Section 2 contains our percolation criteria for the 
existence of infinite avalanches leading naturally to the definition of the functional 3. In Section 3 
we show that 3 C = | is the unique critical "point" of our model, thus ruling out the possibility of 
an intermediate phase. Section 4 proves the above relations for the critical exponents. Finally, in 
Section 5 we develop a coupling argument which is the core of the proofs of the aforementioned 
results in Sections 3 and 4. The principal results of this paper are Theorem 2.4 (Section 2.2), 
Theorem 3.1 (Section 3.1) and Theorem 4.1 (Section 4.1). 

2. Percolation criteria 

2.1 Simple percolation bounds. 

We start by deriving criteria for the presence and absence of an infinite avalanche based on a 
comparison to site percolation on Let x* denote the maximum of the support of p, i.e., 

x* = sup{yG [0,1]: p([y,l}) > 0}, (2.1) 

and let us define 6*5 by 




(2.2) 



It is noted that if X + v < 6^, then the largest value that X a (t) for any a G could conceivably 
achieve (just prior to its own avalanche) is fy,. 

The following is based on straightforward percolation arguments: 

Proposition 2.1 (1) //p([^p, 1]) > f, then ¥ p (\^ \ = 00) > for all v > 1 - x*. 

(2) If either 6 b < 1 or 6 h > 1 and p([l - \6 b , 1]) < \, then W p (\s/W \ = 00) = Ofor all v > 0. 
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In both cases we note that the quantity | on the right-hand side of the inequalities is the 
percolation threshold for T b . Obviously, this is no coincidence; indeed, the proof of part (1) is 
easily generalizable to any transitive infinite graph. 

Proof of Proposition 2.1. Let us start with (1): A site a / is called occupied if X a > 
1 — -j, while the root is called occupied if X + v > 1. Denoting by the connected 
component of occupied sites containing the origin, it is not hard to see that^M D <*fH Indeed, 
assuming X + v > 1, each daughter site of the origin receives at least |; those daughter 
sites er with X CT > 1 — | will be triggered, which will in turn cause avalanches in the next 
generation of occupied sites, etc. Evidently, whenever the occupied sites percolate, there is an 
infinite avalanche. 

Part (2) is proved in a similar fashion. Suppose first that 9 b > 1 and call a site a / occupied 
if X a > 1 — \9 b , and vacant otherwise. The definition is as before for the root. As observed 
previously, if X + v < 9 b , then no site receives more than ^9 b from its parent. Under these 
circumstances, a vacant site will never avalanche and, denoting the occupied cluster of the origin 
by we have st$) C tfM. Since p([l- \9 h , 1]) < £ was assumed, we have that \^\ < oo 
almost surely and thus |^(")| < oo whenever X + v < 9 b . It is then easy to show, however, 
that \&/( v ) | < oo almost surely for all v > 0. Indeed, let k > be an integer so large that 

{ Xir + v- e b )b- k <e b -i. (2.3) 

If a is a site with |c| = k that has been reached by an avalanche, then a could not receive 
more than 

x^b- 1 + • • • + b-( k -V) + b- k (x + V ) = b~ x e h + b- k (x + v-e b ) (2.4) 

from its parent. Now, if a is vacant, then the maximal possible value for X a (k) (i.e., prior to 
its own avalanche) is no larger than 1 + b~ k (x+ + v — 9b)- By (2.3), this amount is strictly less 
than 9 b , so by our previous reasoning, a cannot give rise to an infinite avalanche. By absence 
of percolation, there is a "barrier" of vacant sites above the (k + l)-st layer in T b , that every 
path from the root to infinity must pass through. Our previous arguments show that the avalanche 
cannot go beyond the union of occupied connected components rooted at Hence, \s^^> \ < oo 
with probability one. 

The case 9 b < 1 is handled analogously. Indeed, a simple calculation reveals that the right- 
hand side of (2.4) plus x+ is eventually strictly less than one and the avalanche terminates within 
a deterministic (^-dependent) amount of time. □ 

The arguments in the proof immediately give us the following corollary: 

Corollary 2.2 If9 b ^ 1 and p([l — ^9 b , ^p)) = 0, then there is an infinite avalanche ifandonly 
if occupied sites, i.e., sites a G T b with value X a > 1 — p percolate. In addition, if X + v < 9 b , 
then &f( v ) coincides exactly with the occupied connected component of the root. 

Remark 1. The exceptional cases, 9 b = 1, can only arise from the circumstance that = 1 — |. 
(Notice that the proof of Proposition 2.1(2) does not apply because the inequality (2.3) cannot 
be satisfied.) For 9 b = 1, the situation is marginal and, in fact, slightly subtle. Indeed, if = 



A MODEL OF ORGANIZED CRITICALITY 



7 



1 — ^ and F(X > x+) = \, then the existence of an infinite avalanche depends on the detailed 
asymptotic of ¥(X > x+ — e) as e j 0, see Remark 2 in the next section. We exclude the 
cases 9b = 1 from our analysis because we believe that this "pathological" behavior is in no 
way generic. 



2.2 Phase transition. 

As is seen from Corollary 2.2, in certain cases the avalanche problem reduces to the usual (in- 
dependent) percolation model. The general problem can also be presented as a percolation phe- 
nomenon albeit with correlations. Indeed, let X\, . . . , X n are i.i.d. with distribution p and let 

In the case n = 0, we let Q { ° ] = 9. Similarly, for a G T 6 , we define Q { a 9) by (2.5) with n = \a\ 
and Xi, . . . , X\ a i being the values along the unique path connecting a to the root. Explicitly, we 

set = 9 and define 

QP=x* + \q%)> ° + (2 " 6) 

Note that here 9 plays the role of the quantity X + v. Clearly, Qn = Q a , whenever n = \a\. 

Proposition 2.3 Let v > and let 9 = X + v. For each a 6 T;,, let us call a open if Q a > 1 
and closed otherwise. Then a G &/( v ) if and only if a belongs to the open cluster containing the 
root. In particular, percolation of open sites is the necessary and sufficient condition for infinite 
avalanches. 

Proof. By definition, Q ' = 9 = X +v. Now, if X a {t) = Q a for a site a G Tj that avalanches 
at time t = \a\, then any daughter site a' of a will have its value updated to 

X a ,(t + 1) = X a , + l -Q^ = Q^. (2.7) 

Hence, if the site a G avalanches at time t = \a\, then Q a ' = X a (t) > 1. It follows 
that g/( v \ with v = 9 — X , is the set of sites that are open and connected to the root by a path 
of open sites. □ 

Remark 2. Let us indicate what makes the case x* = 1 — | so subtle. Given a sequence (ck) of 
positive numbers, let us call a G open if X a > x+ — C| cr |6~l cr l and closed otherwise. Letting 
Pk = F(X > x+ — Ckb~ k ) and supposing, e.g., bpu = 1 + k~ 1 / 2 , a general result of Lyons [9] 
implies that the open sites percolate. An easy argument shows that if a is connected to by a path 
of open sites, then C$P > 1 + b~ k (v — 1 — J2e<k c ?) for v = 9 — X . Thus, if v > 1 + ^ fc>0 c&, 
then, by Proposition 2.3, there is an infinite avalanche with a non-zero probability. 

On a similar basis, we can write down the necessary and sufficient conditions for divergence 
of the expected size of avalanches. The criterion will be based on the asymptotic growth of the 
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quantity 

Z n (9) =F(Q[ e) >1, k = 0,...,n), n>0. (2.8) 
Notice that Z n {6) = whenever 9 < 1. 

Theorem 2.4 (1) For all 9 > 1, the limit 

3=3(p)= lim Z n (9) l l n (2.9) 

ra— +00 

exirfj a«(i is independent of 9. 

(2) For all p, p' £ M., the function a i— ► 3(ap + (1 — a)p') is continuous in a £ [0, 1]. 

(3) Let p £ A4 and let x* correspond to p via (2.1). Define } c = 4. If $(p) < 3c> 
Epd^^l) < ooforallv G (0,oo), while if ^(p) > } c , thenE p (\s/^\) = ooforallv > 

Theorem 2.4 defines a free-energy like functional 3 and gives the characterization of the diver- 
gence of x^ v \ as already discussed in Section 1.3. The continuity statement in part (2) indicates 
that the sets of "avalanching" and "non-avalanching" measures p G M are separated by a "sur- 
face" (i.e., set of codimension one) of phase transitions. We will not try to make the latter more 
precise; our main reason for including part (2) is to have an interpretation of the limit 3 (p) — ► 3 C , 
which will be needed in the discussion of the critical behavior. Under additional mild restrictions 
on p, it will be shown in Section 4 that E p (|=e/^|) = 00 even for the critical measures p, i.e., 
those satisfying 3(p) = 3 C . 

Proof of Theorem 2.4(1). We will start with the cases 6 = 1 and > Ob which are amenable to 
subadditive-type arguments. Examining Z n+m (0), we may write (by conditioning on Xi, . . . , X m ) 

m 

Z n+m (8) = E(z n (QW) II 1 W f )>i})- < 2 - 10 ) 

3=0 3 

Since 9 ^ Qn^ is manifestly non-decreasing in 6, so is the event on the right-hand side of (2.8) 
and also Z n (9) itself. Notice that if 9 > 9 b , then < 9 for any k > 0, while if 9 = 1, then the 
conditions in (2.8) force Q k > 1. Thus, for 9 = 1 we obtain the submultiplicative bound 

(1) > Z n (l)Z m (l), (2.11) 
while for any 9 > 9b we get the supermultiplicative bound 

(9) < Z n (9)Z m {9). (2.12) 

By standard theorems, Z n (V) x l n tends to a limit, 31, while Z n (9) l / n for 9 > 9b tends to a (pos- 
sibly ^-dependent) limit } e . Moreover, (2.10) in fact implies that Z n+m (9) < Z n (9 b + 9b~ m ) 
and 3e is thus constant for all 9 > 9\>. We will use 3* to denote the common value of $0 for 9 > 9 b . 
Note that ^„(l) 1/n < 31 while Z n {9f/ n > 3^ for all n > 1 and all 9 > 9 b . 

Since 9 Z n {9) is non-decreasing, to prove (2.9), we just need to show that 3+ equals 31. 
If < 1 — |, then 3* = and there is nothing to prove, so let us suppose that > 1 — | for the 
rest of the proof. As it turns out, we will have to address a number of distinct cases. These are 
determined by whether the inequality in x± > 1 — | is strict or not and by whether the quantity 

n e = p([x* - e,x+\) (2.13) 
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is strictly less than $ 1 or not for some (particular) e > 0. Specifically, we will distinguish the 
following cases: 

CASE 1: x* > 1 — £ and k 6 < 31 for some e > with x* — e > 1 — ^. 
CASE 2: x* > 1 - | but K e = 31 for all e > with x* - e > 1 - \. 
CASE 3: x* = 1 - \. 

As is easily observed, CASE 2 represents the situation where p assigns no mass to the interval 
(1 — l,x+), while CASE 3 corresponds to the similar situation when this interval itself is empty. 
In view of the trivial inequality 31 > p([l — ^, x*]), we must eventually have k € < 31 whenever p 
has any mass in (1 — ^, x*). Hence, the first situation is clearly generic. 

In order to address the first two cases (with x* > 1 — \) we need to establish an inequality 
between Z n (9) and 3™ for all 9 £ [1, 9b)- Explicitly, we claim that for x* > 1 — ^ and any 9 G 
[1, 9b), there is an H(9) < 00 such that 



Z n {9) < H(0)fi, 
Indeed, let e > be such that 9b — 9 > e-^ and x* 



+ 



n > 1. 

- e > 1 — ^ and pick m so that 
1 1 1 



,m-l 



H > 

6 m ~ 



(2.14) 



(2.15) 



Consider the formula (2.8) for Z n+m (l) but with the first m coordinates restricted to the event £ = 
{Xi, . . . , X m > x+ — e}. Notice that on £, the conditions involving . . . , are automat- 
ically satisfied. By a derivation similar to (2.10-2.11) we have 

Z n+m {l) > CZ n (9). (2.16) 

Along with the upper bound Z n+m (l) < 3™ +m , this implies (2.14) with H{9) = {li/n t ) m . (We 
note that, since x* is the suppremum of the support of p, we have n e > for all e > 0.) 

Now we are ready to prove that 3* = 31 in all of the three cases above: 
CASE 1 : Suppose that x* > 1 — | and n e < 31 for some e > with x* — e > 1 — Let 1 
be small enough that 9 e = x* — e + ^ < 9b- Then 

Z n -i(9 e ) 



> 



Z n {9) < K e Z n ^{9) + (1 - K e )Z n ^(9 t ) = Z n _i(0) 
Using (2.14) and the bound Z n _i(0) > 3™" 1 we obtain 



Z„(0) < Z n ^(9) 



K e + (1 - K £ )ff( 



« e + (l 



(*)J 



n-1 



(2.17) 



(2.18) 



Let K e (n) denote the quantity in the square brackets, and let us set n = 2m in (2.18) and iterate 
the bound m times. This gives Z2 m {9) < K e (m) m Z m (9). If we still entertain the possibility 
that 31 < 3*, then the m — > 00 limit gives 3* < linim^oo K e (m) = K e , which contradicts the 
bound 31 > K e . Therefore, once x* > 1 — \ and k € < 31 for some e > we have 31 = 3*. 

CASE 2: Suppose now that x* > 1 — \ but k 6 = 31 for all e > with x* — e > 1 — ^. 
Notice that this in fact implies that 31 = p({x^}). We first observe, using (2.16) with n = 1, 
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that Z m+ i(l) > n™p([l — j;6,x±\) whenever 9, e and m satisfy (2.15). As a consequence of 
(2.11), we have 



3i > 



P ([i-le,x*]) 



i 



m + l 

(2.19) 



Now if p((l — l9 b ,x+)) > 0, we would have p((l — l9 b ,x±\) > = K e which would 

by (2.19) imply that fa > re e , a contradiction. (This fact will be important later, so we restate 
it as a corollary right after this proof.) Hence, we must have p((l — ^06, x*)) = 0. To prove 
that 31 = 3*, let 6> > b be small enough that 9q = I + \{9 — 9 b ) < 9 b . Now either X fc = 

for all fc = 1, . . . , n, or there is a A; such that < 1 — ^9 b . Noting that then < 9q, we 
thus have 



Z n {9) < p{{x*}) n + J2p{{x*}) k V([0,1- £0 6 ])Z n _ fc (0 o ). ( 2 - 20 ) 
fc=i 



Using (2.14), this gives Z n {9) < + 713™ V([0, 1 - \9 b \)H(9o), proving that 3* = 31 holds 
for CASE 2 as well. 

CASE 3: Suppose finally that = 1 — £ and note that then 0& = 1. Immediately, we 
have Z n (l) = p({x+}) n and therefore 31 = p({x*}), while for any 9 > 9 b we have n^oiQ/if ^ — 
1} C fTk=i{X k > - b~ k {9 - 9 b )}. Therefore, 

n 

Z n {9)< l[F(X k >x*-b- k (9-9 b )), (2.21) 

k=l 

which implies that 3* < lim^^ F{Xi > x+ - b~ k (9 - 9 b )) = p({x*}) = 31. □ 

This completes the proof of part (1) of Theorem 2.4. As mentioned earlier, we would like to 
underscore one aspect of the above proof. 

Corollary 2.5 Let p G M. and suppose that x± > 1 — | and p((l — ^9 b , x*)) > 0. Then there 
is an e > with x* — e > such that ^(p) > /o([x* — e, x*]). 

Proof. See the argument following (2.19). □ 

Next we will prove the continuity of a 1— ► 3(ap + (1 — a)p') as stated in Theorem 2.4(2): 

Proof of Theorem 2.4(2). Throughout this proof we will write Z^(6) instead of just Z n (0) 
to emphasize the dependence on the underlying measure p. Let po,pi € M. and let p a = 
(1 — a)po + ap\. Clearly, to prove (2), it suffices to show that a ^ l{p a ) is right continuous 
at a = 0. 

Fix a > and let (T k ) be a sequence of 0, 1-valued i.i.d. random variables with Prob(Tfc = 
0) = a. Let (X k ) and (X' k ) be two independent sequences of i.i.d. random variables, both 

independent of (T k ), with distributions p$ and p^, respectively. Let (X^) be the sequence 
defined by 

X ( h a) = T k X k + (1 - T k )X' k , k > 1. (2.22) 



A MODEL OF ORGANIZED CRITICALITY 



11 



Clearly, (X™) are i.i.d. with joint distribution p^. Let us use P a to denote the joint distribution 
of(X k ),(X' k ),aad(T k ). 

Let Q^f' a) be given by (2.5) with X 1 ,...,X n replaced by x[ a) , X { n ] . Then zi Pa ] (9) is 
given by (2.8) with Qn^ replaced by Qn' a ^ and P replaced by F a . As will be seen shortly, the 
main object of interest is the conditional expectation given the values (T k ): 

Z n , a {6\(T k )) =F a {Qf' a) > l,£ = 0,...,n\(T k )). (2.23) 

Indeed, let 9 G [1, 9b] and, given (Tk), let (ij) be the connected blocks of sites A; G {0, . . . , 1} 
such that Tfc = 1 and let (Jj) be the connected sets of sites with T k = 0. By (2.22), the xj^ 
for k G Jj are distributed according to po, while those for A; G Jj are distributed according to pi. 
Then an analogue of (2.10) for the quantity in (2.23) along with the bounds Z n (l) < Z n (9) < 
Z n {9b) f° r £ [lj 6b] allow us to conclude that 

i 3 i 3 

In order to estimate the right hand side of (2.24), note that the existence of the limit in (2.9) 
implies that for all S > there is C$ G [1, oo), such that for both p = po and p = pi, 

Z ( n p \9 b ) < C s (l + ^(pf, n > 1- (2-25) 

Let E a denote the expectation with respect to F a . Using (2.25) in (2.24), observing that the total 
number of occurrences of Cs is less than 2k\{Y), where k\{Y) = J2j \Jj\> an d noting that ki(Y) 
has the binomial distribution with parameter a under P a allows us to write 

Z^(9)=E a Z n>a (9\(T k )) < (l + 6) n ((l-a)}( Po ) + aCfo( Pl )) n . (2.26) 

By taking n — ► oo, we get \\m. a ^l(p a ) < (l+<5)3(Po)- But (5 was arbitrary, hence, lim^o 3 < 
3(po). The argument for the lower bound, lim^o 3(p a ) > 3(po), is completely analogous. □ 

Finally, we also need to prove part (3) of Theorem 2.4: 
Proof of Theorem 2.4(3). By Proposition 2.3, a G st^ is exactly the event that the path between 
(and including) a and consists of sites a' with Q^) > 1, where 9 = X + v. But then 

F p (a G s/W) = E p (Z ]a] (X + v)), (2.27) 

where the final average is over X . To get the expected size of srf( v \ we sum over all a, 

MW {V) \) = YsM Z n{X + v))b n . (2.28) 

n>0 

The existence of the limit Z X J n (9) independent of 9 (for 9 > 1) tells us that Ep(\*/W\) < oo 
whenever 3(p) < 3 C , while E p (|i2/^|) = oo once l(p) > } c and v > 1 — x±. □ 



3. Absence of intermediate phase 



12 



M. BISKUP, PH. BLANCHARD, L. CHAYES, D. GANDOLFO AND T. KRUGER 



3.1 Sharpness of phase transition. 

The goal of this section is to show that the phase transitions defined by presence/absence of an 
infinite avalanche and divergence of avalanche size occur at the same "point," % c = \. This rules 
out the possibility of an intermediate phase. Moreover, we will prove that the transition is second 
order in the sense that there is no infinite avalanche at 3 = 3c- 

Unfortunately, our proof will require certain restrictions on the underlying measure p. The 
delicate portion of [0, 1] is the region I = [1 — ^6 b , ^-). Clearly, some conditions are needed to 
ensure that there is not too much mass at the right-end of I — i.e., that p([(l — \ — e, 1 — -|)) —>■ 
sufficiently fast as e [ — to avoid the sort of counterexamples described in Remarks 1 and 2 of 
Section 2. With a lot of additional work than what is about to hit, all of the forthcoming can be 
proved under the assumption that p has an L p density, for a p > 1, in the interval /. However, 
this requires dealing with "singularities" in the region above /. (The region below / is of no 
consequence because any directed path in the avalanched set can only harbor a finite number of 
values from this set.) Notwithstanding, most of the interesting mathematics — with only a fraction 
of unpleasant technicalities — is captured by assuming that the measure p has a bounded density. 

Definition. Let M}' be the set of Borel probability measures p on [0, 1] that are absolutely 
continuous with respect to Lebesgue measure on [0, 1] with the associated density (j) p bounded 
in L°° norm on [0, 1], i.e., \\4> P \\oo < °°» an d that obey p{[l — |, 1]) > 0. 

The requirement that p has no positive mass in [1 — |, 1] represents no additional loss of 
generality since the opposite case, namely x* < 1 — \, actually has 3(p) =0 and is therefore 
far away from having an avalanche (see Theorem 2.4). It is worth noting that Ai^ is a convex 
subset of M. The ability to take convex combinations of elements of will be crucial in the 
discussion of the critical behavior, see Section 4. 

Our second main theorem is then as follows: 

Theorem 3.1 Suppose that p G and define 3 C = ^. 

(1) If lip) < 3c, thenP p (\£/^\ = 00) = for all v e (0, 00). 

(2) Ifl(p) > 3c, thenP p (\£/^ \ = 00) > for all v G (l-x*,oo). 

The proof of Theorem 3.1 requires introducing two auxiliary random variables V^o and Qoo. 
These will be defined in next two subsections, the proof is therefore deferred to Section 3.4. The 
random variable Qoo will be a cornerstone of our analysis of the critical process, see Section 4. 
The underlying significance of both Voo and Qoo is the distributional identity that each of them 
satisfies. 

3.2 Definition of V^. 

In this section we define a random variable Voo which is, roughly speaking, the minimal value 
of v that needs to be added to the root in order to trigger an infinite avalanche. For n > 1 let 



V n = \ni{v e (0,oo): / X( v \t), t = 0,...,n- l}. 



(3.1) 
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(A logical extension of this definition to n = is Vq = 0.) In plain words, if v > V n , then the 
avalanche process will propagate to at least the n-th level. Clearly, V n is an increasing sequence; 
we let Voo denote the n — ► oo limit of V n . Formally, Voo can be infinite; in fact, since the 
event {V^ < oo} is clearly a tail event, P p (V r 00 < oo) is either one or zero. 
Let us use ^ n to denote the distribution function of V n , i.e., 

tf n (0) =F p (V n <0). (3.2) 

The aforementioned properties of V n lead us to a few immediate observations about First, ^> n 
is a decreasing sequence of non-decreasing functions. Second, the limit 

= lim *„(<?), (3.3) 

n— >oo 

exists for all ■& G (0, oo) and = P p (V r 00 < #). Third, we have * ^ if and only 

if Pp(V 00 < oo) = 1. Moreover, each of Vl/ n is in principle computable: 

Lemma 3.2 Let p £ M. Then the sequence (^ n ) satisfies the recurrence equation 

* n+1 ( 1 ?)=E^$ 6 (* n (^))l {X0 > 1 _^, n>0, (3.4) 
where ^o(^) = "l{tf>o} an d 

<S> b (y) = l-(l-y) b , 0<y<l. (3.5) 

Proo/ Let denote the subtree of ¥5 rooted at a and let denote the random variable 
defined in the same way as V n but here for the tree . Then we have 

{Vn+i < = {X > 1 - tf} n { min < ^r^j. (3.6) 

I ae{i,...,b} b ) 

But for all a € {1, ... , 6}, the 's are i.i.d. with common distribution function * n , so we have 

P„( min V n W < 0) =$ 6 (¥„(0)). (3.7) 

(7G{l,...,fe} 

From here the claim follows by noting that are independent of X . □ 
Corollary 3.3 Le? p £ M.. Then the distribution function of Voo satisfies the equation 

= Ep($6(^(^))l { x >i-^}), ■d > 0. (3.8) 

/Voo/ This is an easy consequence of (3.4) and the Bounded Convergence Theorem. □ 

On the basis of (3.8) and some percolation arguments, the answer to the important question 
whether \t = or not can be given by checking whether ^(1?) = for reasonable values of 

Proposition 3.4 Let p <E M^. Suppose that ^ ^ 0. Then 

inf{tf> 0: > 0} = (3.9) 
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Proof. Let denote the infimum on the left-hand side of (3.9). Note that x* > 1 — \ by p G M^. 
Since p is absolutely continuous with respect to the Lebesgue measure on [0, 1], there is an 77 > 
such that x* — 77 > 1 — | and /?([x* — r/, x*]) < 4. Now | is the threshold for the site percolation 
on Tb, so the sites with X a > x* — r\ do not percolate. Let G„ = {cr G T^, : |<r| = n} be the n-th 
generation of T&. Pick two integers A/", AT' such that N' > N and let Hn,N' be the event that G^v 
and Ga?/ are separated by a "barrier" of sites a with X CT < x+ — r). By taking N f S> AT s> 1, the 
probability of Hn,n' can be made as close to one as desired. 

Let ■!? > and pick A?o so large that #6"^° is less than \. Find N, N' > N such that 
1 - P P (H N ,N') is strictly smaller than P p (\*/W\ = oo), i.e., we have P p ({|.c/W| = oo} n 
7~Cn,N') > 0. Now for any e G (0, we will produce a configuration with an infinite avalanche 
that has a starting value v = 1 — x* + e. Draw a configuration subject to the constraint 
thatX CT > x*-e for allcr G T b with|cr| < A/"'. Let belong to the set {\^^\ = oo}r\H N ,N> 
and define X' a by putting 

K = \x^x„ if WS iV, (3]0) 

IX CT , otherwise. 

Let A^ (t,) (t) denote the process corresponding to the initial configuration (AT^) and initial value 
v > 0, and let (t) be the corresponding process for (X a ) and Let ^''^ and be the 
corresponding avalanche sets. 

The configuration (X a ) exhibits an infinite avalanche, so there is a site a on one of the afore- 
mentioned "barriers" separating Gat and Gn 1 , which belongs to an infinite oriented path in- 
side g/^. By the assumption that x* — 7] > 1 — ^ it is clear that, if v > 1 — x* + e and t = \a\, 
then srf'^ will reach a. But X' a , > X a for all sites on the path from to a, so we have 

X'M(t) -lf)(t)>,- £ -^> 0, (3.11) 

where we used that b N $ < \ and e < ^ to derive the last inequality. Now the set si^ contains 
a path from a to infinity and, by (3.11) and X' a , > X a / for a' "beyond" a, this path will also 
be contained in s^''^ v \ Consequently, an infinite avalanche will occur in configuration (X' a ) 
starting from a value v > 1 — x* + e whenever it did in configuration {X a ) starting from This 
establishes = 1 — x*, as claimed. □ 



3.3 Definition of Qoo. 

The second random variable, denoted by is a limiting version of the objects Q n defined 
in (2.5). Let Y = (Yi, Y2, . . . ) be a sequence of i.i.d. random variables with joint distribution 
P = p N . These are, in a certain sense, the same quantities as the X's discussed earlier, however, 
the F's will be ordered in the opposite way. Similarly to (2.5), let 

<i = ^ + ^W--- + ^n + ^ TT , l<k<n. (3.12) 

For completeness, we also let \ = 9. 
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Let B be the Borel cr-algebra on [0, 1] N equipped with the standard product topology. Suppose 

that p([^j±, 1]) > 0— which is assured if p € M b . For any n > 1 and > 1, let be the 
conditional law on B defined by 

Pf(.) =P (.|QW> M = 2 ,...,n), (3.13) 

The latter is well defined because P(Qj$ > 1) > for all £ = 2, . . . , n, {qJJ ^ x i 316 
increasing and P(-) is FKG. Intentionally, the variable Y\ is not constrained by the conditioning 
in (3.13). 

Next we give conditions for the existence of the limiting law lim^oo P„ : 

Proposition 3.5 Let p G M b and let 9q > 9b. Then there exists numbers A = A(p, 9q) < oo 
and ( = ((p) > such that for all bounded measurable functions f = /(Yi, . . . , Y^) and all 
0,0'e[l,0 o ], 

|eW i(/ ) -E(?')(/)| < Ae-^WflU n>k. (3.14) 
In particular, whenever 9 > 1, Z/m/f /aw 



n^oo 



P(-) = lim (3.15) 



exists and is independent of 9. Moreover, the quantities A(p,9q) and £(p) are bounded away 
from infinity and zero uniformly in any convex set M C M b with finitely many extreme points. 

The proof of Proposition 3.5 uses a coupling argument, which requires some rather exten- 
sive preparations and is therefore deferred to Section 5. (The actual proof appears at the end of 
Section 5.3.) 

We will use E to denote the expectation with respect to P whenever the latter is well defined. 
Let us define a random variable Qoo on ([0, 1] N , B, P) by the formula 

k>i 

Notice that Qoo is supported in [|, 9j,], because Y\ is not constrained by the conditioning in (3.13). 

Corollary 3.6 Let p G M b and let 9 > 1. Let \ be as in (3.12), where the variables 

Yi, . . . , Y n are distributed according to Pn \ Then \ tends to in distribution as n — > oo. 
Moreover, for each 9q > 9b and each C < oo there are constants D = D(p, 9q) < oo and c = 
c(p) > such that if f(9) is a function obeying the Lipschitz bound on [0, 9q], 

\f(9)-f(9')\ <C||/|U0-0'|, 0,0'e[O,0 o ], (3.17) 

where ||/||oo = sup e < e(J \f{9)\, then 

(/(<})) -E(/(Qoo))| < H/llooe-^ (3.18) 

holds for all 9 £ [1, 9q\. The quantities D(p, 9q) and c(p) are bounded away from infinity and 
zero uniformly in any convex set M C M b with finitely many extreme points. 
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The proof of Corollary 3.6 is given in Section 5.4. As already mentioned, a principal tool for 
our later investigations will be the distributional identity for Qoo stated below. 

Proposition 3.7 Let p G M\ If X is a random variable with law P = p, independent 
ofQoo, then 



Qoo>l =P(Qoo€-)- (3.19) 



The proof of Proposition 3.7 will also be given in Section 5. Proposition 3.5 and the proof 
of Proposition 3.7 immediately yield an extension of Theorem 2.4(1), stated as Corollary 3.8, 
which will also be useful in subsequent developments. The proof of Corollary 3.8 is given in 
Section 5.4. 

Corollary 3.8 Suppose that p G M^. Then %{p) = P(Qoo > 1)- Moreover, the limit 

i> p {9) = lim Z n (e)l{p)- n (3.20) 

exists for all 9 > and, for all 9q > 9b, there are A' = A'(p, 9q) < oo and Q 1 = C'(p) > 
such that 

\Z n {6)l(pY n - MV)\ < A 'e~ Cn (3-2D 
holds for all 8 G [0, #o] and all n > 1. Furthermore, the function ip p has the following properties: 

(1) WW e (0, oo) for all 9 > 1 while ip p {9) = Ofor 9 < 1. 

(2) 9 i— ¥ ip p {9) is non-decreasing and Lipschitz continuous for all 9 > 1. More precisely, 
there is aC = C(p,9 ) < oo such that \i/j p {9)-tp p {9')\ < Cip p {9 )\9 -9'\for all 9,9' G 
[l,<9o]. 

(3) If p,p' G and p a = (1 — a)p + ap' for each a G [0, 1], then a i— > ip Pa (9) is 
continuous in a G [0, l]/or a// # > 0. 

77ie quantities A'(p, 9q), C'(p) an d C(p, 9q) are bounded away from infinity and zero uniformly 
in any convex set M C Ai^ with finitely many extreme points. 

Remark 3. The Lipschitz continuity of 9 i— > V'p(^) i s a direct consequence of our assumption 
that p has a bounded density 4> p with respect to the Lebesgue measure on [0, 1]. If 4> p is only 
in L p ([0, 1]) for some p > 1, then the appropriate concept will be Holder continuity with a p- 
dependent exponent. The same will be true for various other Lipschitz continuous quantities later 
in this paper. 

3.4 Proof of Theorem 3.1. 

With random variable Qoo at our disposal, the sharpness of the phase transition in our avalanche 
model is almost immediate. 

Proof of Theorem 3.1. Let p G and abbreviate 3 = l{p). Let x* be as in (2.1). We begin by 
introducing the quantity 

G n = E(^%)l {Qoo > 1} ). (3.22) 
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The recursive equation (3.4) and Proposition 3.7 then give 
= 3E($Jf n (^)1 {QK > l} )) l 



where we have used the fact that 3 = P(Qoo > 1) from Corollary 3.8. 

Let us first analyze the cases 63 < 1. By using Jensen's inequality in (3.23) we get that 

G n+1 < 3 MGn) < iMGn). (3.24) 

An inspection of the graph of y Q^y) reveals that if (3.24) holds, then G n — > 0. By Proposi- 
tion 3.4, this is compatible with ^ ^ only if 

Qc 



;oo 



< 1 - x* P-almost surely. (3.25) 



However, a simple argument shows that esssupQoo = x *j^i whenever 3 > 0. This contradicts 
(3.25), because > (as implied by p G M. b ) forces 1 — x+ < x^-^j. Thus, if 63 < 1, 
then must be identically zero. 

Next we will attend to the cases 63 > 1. We will suppose that ^ n — > and work to derive 
a contradiction. Since n h-> ty n is a monotone sequence of monotone functions, the conver- 
gence to ^ is uniform on [0, 1] and, in particular, on the range of values that ^Qoo takes. Using 
that $b(y) >by — \b{b — l)y 2 for all y <E [0, 1] and invoking (3.23), we can write 

G n+l >&3(l-e n )G n , (3.26) 

where e n = \{b — l)\I/ n (l). Since 63 > 1 and e n — > 0, we have G n+ i > G n for n large 
enough. An inspection of (3.4) shows that, since x* > 1 — \, we have \I/ n (i?) > for all 1? > \. 
Hence G n > for all n > 0. But then (3.26) forces G n to stay uniformly bounded away from 
zero, in contradiction with our assumption that G n — > 0. Therefore, once 63 > 1, we must 
have * 0. □ 



4. Critical behavior 

4.1 Critical exponents. 

In this section we establish, under certain conditions on p, the essential behavior of the model at 
the critical point 3 C = ^. In particular, we describe the asymptotics for the critical distribution of 
avalanche sizes, the power law behavior for the probability of an infinite avalanche as 3 J. 3 C and, 
finally, the exponent for the divergence of as 3 t 3c- 

Theorem 4.1 Let p G M b and let x± be as in (2.1). Suppose i(p) = 3 C , where 3 C = i. Then 
there are functions t : T : (1 — x*, 00) — > (0, 00) and : [0, 00) — > [0, 00) swc/i ?/zatf the following 
holds: 
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(1) If p' G Ai^ and p a = ap' + (1 — a)p satisfies l(p a ) < J c for all a G (0, 1], then for 
all v > 1 — x*, 

E ^(l^ (w) l) = 3c -3(p a ) [ 1 + o(1) J' 

(2) For all v > 0, 

P p (|^ (,;) |>n) = ^[1 + (1)], n^oo, (4.2) 
where 0(v) > 0/or v > 1 — x*. 

(5) If p' G .M a«(i jO a = ap' + (1 — a)p satisfies l{p a ) > far all a G (0, 1], then for 
all v > 1 — x*, 

PpJI^I =oo) = ( 3 (p a ) - 3c ) T( V ) [1 + (1)], ajO. (4.3) 

Remark 4. The proof of Theorem 4. 1 makes frequent use of the properties of the random vari- 
able Qoo defined in Section 3.2. The relevant statements are Propositions 3.5 and 3.7 and Corol- 
laries 3.6 and 3.8, whose proofs come only in Section 5. Modulo these claims, Section 4 is 
essentially self-contained and can be read without a reference to Section 5. 

Part (1) of Theorem 4.1 can be proved based on the already-available information; the other 
parts will require some preparations and their proofs are postponed to the next section. 

Proof of Theorem 4.1(1). Letp, p' G be such that l(p a ) <% c =l(p)fovp a = (l — a)p + ap' 
and all a G (0, 1]. Let x (w) («) = e pA\^ {v) \)- b Y ( 2 - 2 8), 

X {v) {a) =^E Pa (4^)(X + *,))&", (4.4) 

n>0 

where E Pa is the expectation with respect to X in p a and Z?f ^ is defined by (2.8) using p a . 

In order to estimate the sum we will use A" and (" to denote the worst case scenarios for the 
quantities A'(p a , 9q) and C'(Pa) from Corollary 3.8. Explicitly, we let A" = sup < a <i 
and (" = info< Q <i C(p a ), where 6q > @b is to be determined shortly. Note that A" < oo and 
C" > by uniformity of the bounds on A'(p a , 6$) and ('(p a ) in the convex set N = {p a : a G 
[0, 1]}. Then we have, for all n > 1 and all 6 G [1, O ], 

6"i#»>(0) = 6"3(Pa) n ^ a (^) + 6«j(p a )"i5; n (0), (4.5) 

where i/) Pa (0) is as in (3.20) while E n (8) is the "error term." Using the bounds from Corol- 
lary 3.8, E n (0) is estimated by \E n (9)\ < A"e~^" n . By continuity of a i-> ip Pa (6), we get 

EWe) = Mij^, (4 . 6) 

n>0 1 

where o(l) tends to zero as a j uniformly on compact sets of 9 G [1, Oq]. 

Let t(v) = b~ 1 ~E p (ip p (X0 + v)) and note that r(u) > for all v > 1 — x*. Let us take the 
maximum of x* + v and 26*5 for the quantity 9$ above. Then (4.4) and (4.6) imply 
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where o(l) tends to zero as a j 0, for all v > 1 — x*. □ 

It remains to establish parts (2) and (3) of Theorem 4.1. To ease derivations, instead of looking 
at the asymptotic size of &/( v \ we will focus on a slightly different set: 

>}, if^(^) = 0, 

(^{o- G T 6 : m(cr) G js^"^}, otherwise. 

(Here we take a/^'^ = whenever 9' < 1.) Clearly, is the original avalanche set to- 
gether with its boundary (i.e., the set of sites in where the avalanche has "spilled" some 
material). Since both sets are connected and both contain the root (with the exception of the 
case £/( 9 ~ X0 ^ = 0), their sizes satisfy the relation: 

\m {e) \ = {b-l)\^ {6 - X0) \ + l. (4.9) 

(This relation holds even if ,e/( e ~ X0 ) = 0.) The asymptotic probability of the events {J^") | > 
n} as n — > oo is thus basically equivalent to that of > (b — l)n}. 

4.2 Avalanches in an external field. 

Following a route which is often used in the analysis of critical systems, our proof of Theorem 4. 1 
will be accomplished by the addition of extra degrees of freedom that play the role of an external 
field. Let A G [0, 1] be fixed and let us color each site of Tf, "green" with probability A. Given 
a realization of this process, let denote the random set of "green" sites in T&. Let Po,a(") be 
the joint probability distribution of the "green" sites and (X a ). The principal quantity of interest 
is then 

B 00 (6»,A)=P p>A (^nSf ^0). (4.10) 

It is easy to check that, as A J, 0, the number B oc> (9, A) tends to the probability W p {\B§^ \ = oo). 
In particular, Theorem 3.1 guarantees that B oc (9, A) — ► as A [ if 3 (p) < 3 C , while B^O, A) 
stays uniformly positive as A J, when 3 (p) > } c and 9 > 1. 

Let ip p (0) be as in Corollary 3.8 and let c p G (0, 00) be the quantity defined by 

_L = b -l^ [HMX+ ^))]>\ > l). (4.11) 

Here X and Qoo are independent with distributions P = p and P, respectively. It turns out that 
the asymptotics of B^O, A) for critical p can be described very precisely: 



Proposition 4.2 Let p G M satisfy %{p) = 3 C . For each 9 G (0, 00), 

Bqq{9, A) 

MO 



lim =c P M0)- ( 4 - 12 ) 



Proposition 4.2 is proved in Section 4.4. Now we are ready to prove Theorem 4.1(2): 
Proof of Theorem 4.1(2). We begin by noting the identity 

f?oo(g ' A) = ^P p (|^ (9) | > n)(l - A)"" 1 , A G (0,1], (4.13) 



n>l 



20 



M. BISKUP, PH. BLANCHARD, L. CHAYES, D. GANDOLFO AND T. KRUGER 



which is derived by expressing P p (|^W| = n) as the difference between P p (|^W| = n) and 
Fpd^l = n+1). Since B^O, A) = \f\ {c p ^ p (6)+o(l)) as A | and since n ^ P p (|#W | > 
n) is a decreasing sequence, standard Tauberian theorems (e.g., Karamata's Tauberian Theorem 
for Power Series, see Corollary 1.7.3 in [3]) guarantee that 

P p (|^ (e) | >n) = c p ^jI^= [1 + o(l)] , n ^ oo, (4.14) 
r(i) \Jn 

(Strictly speaking, the above Tauberian theorem applies only when tp p (9) > 0; in the opposite 
case, i.e., when 6 < 1, we have SS 1 ^ 1 = {0} and there is nothing to prove.) In order to obtain the 
corresponding asymptotics for Ppd^/^j > n), we first note that, by (4.9), 

>n) = ¥ p (\^ X0+v ^\ > (b- l)n + l). (4.15) 

By applying (4.14) on the right-hand side and invoking the Bounded Convergence Theorem, we 
immediately get the desired formula (4.2) with 

& ^ = (p-iy^ MMXz + v)), (4.16) 

where E p is the expectation over X . Clearly, v 0(v) is non-decreasing because 9 ^ ipp{9) 
is non-decreasing, while 0{v) > for v > 1 — because tp p (9) > for 9 > 1. □ 

Similarly we can also describe the asymptotics of P p (|^^| = oo) as i(p) [ $ c : 

Proposition 4.3 Let p, p' G and define p a = (1 — a)p + ap'. Suppose that i(p) = } c 
and $(p a ) > 3 C for all a 6 (0,1]. Then for all 9 £ (0, oo), 

p p(l^ ) l = 00 ) = ^ p(g ) + (l) > a i0 , (4.17) 

where ip p (9) is as in Corollary 3.8 and c p is as in (4.11). 

Proposition 4.3 is proved in Section 4.5. Now we are ready finish the proof of Theorem 4.1(3): 
Proof of Theorem 4.1(3). By (4.9) we clearly have that 

P p (|^)| = oo) = F p (\^ X0+v ^\ = oo). (4.18) 

By conditioning on X + v = 9 and invoking (4.17), we can easily derive that the asymptotic 
formula (4.3) holds with T given by T (v) = bc 2 p E p (i> p (X + v)). □ 

As we have seen, Propositions 4.2 and 4.3 have been instrumental in the proof of Theo- 
rem 4.1(2) and (3). The following three sections are devoted to the proofs of the two propo- 
sitions. After some preliminary estimates, which constitute a substantial part of Section 4.3, we 
will proceed to establish the critical asymptotics (Section 4.4). The supercritical cases can then 
be handled along very much the same lines of argument, the necessary changes are listed in 
Section 4.5. 



A MODEL OF ORGANIZED CRITICALITY 



21 



4.3 Preliminaries. 

This section collects some facts about the quantity -Boo(#, A) and its 9 and A dependence. We 
begin by proving a simple identity for B^d, A): 

Lemma 4.4 Let p G M and let $ b be as in (3.5). Then 

B oo (e,\) = \ + (l-\)-\ {e > 1} <!> b (E p B oo (X + le,\)), (4.19) 

Proof. If < 1, then B^^, A) = A and (4.19) clearly holds true. Let us therefore suppose 
that > 1. Let denote the object for the subtree of T 6 rooted at a. Then 

n y + 0} = {0 g ^} u ("{0 ^} n J n sf ^ 0}Y (4.20) 

^ (7=1 ' 

The claim then follows by using the independence of the sets in the large parentheses on the right 
hand side of (4.20) under the measure F Pt \(-). □ 

Our next claim concerns continuity properties of B OQ (9, A) as a function of 9: 

Lemma 4.5 For each p G Ai^ satisfying %(p) < $ c e and each 0$ > 9^ there is a C = 
C(p, 0q) < oo such that 

\B OO (9,X)-B OO (0',X)\ < CB oo (e ,\)\0-0'\ (4-21) 

for all A > and all 9, 9' G [1, 9q\. The bound C(p, 9q) < oo is uniform in any convex set J\f C 
{p G Ai^ : %(p) < } c e} with finitely many extreme points. 

Proof. Let us assume that 9 > 9'. To derive (4.21), we will regard B 00 (6,X) and B 00 {9',X) 
as originating from the same realization of (X a ) and the "green" sites. Then A = B 00 (9, X) — 
Boo(0', X) is dominated by the probability (under F p> \) that there is a site a G T&, a / 0, with 
the properties: 

(1) Q { p > 1 for all a' = m k (a) with k = 1, . . . , |<r|. 

(2) Q ( f ] < 1 but Q ( ° ] > 1. 

(3) ^ o) n # / 0, where is the set for the subtree T ( ° ] rooted at a. 

Indeed, any realization of (X a ) and the "green" sites contributing to A obeys ^ e "> n <S = 
and n ^ 0. But then there must be a site a on the inner boundary of £3^'^ where the 
avalanche corresponding to 9' stops but that corresponding to 9 goes on. (Since 9,9' > 1, we 
must have a / 0.) Consequently, Q y , > 1 for any a on the path connecting a to the root, 

but Qa ' < 1 < Qa ', justifying conditions (1) and (2) above. Since Q K o ' < 0o, and since 
the ^-avalanche continuing on from a must eventually reach a "green" site, we see that also 
condition (3) above must hold. 
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Let p € A4 b be such that j(p) < $ c e. Using the independence of the events described in (1), 
(2) and (3) above, and recalling the definitions (2.8) and (3.13), we can thus estimate 

A < Boo^o, A) Z W\-i@)V\Z\(Qw,i * 1 > Q\Z\,iY ( 4 - 22 ) 

CT eTf,\{0} 

Abbreviate K n (9, 9') = ¥^\q^\ > 1 > Q { °J)- Since Y x is independent of all the other Y\ in 
the measure F n , we have 

K n (6, 6') = {p([l -£,l-£)) :?? '-tf<|0- 9'\b~ n+1 }. (4.23) 

Here resp., play the role of Q^\, resp., Q^2 an d tne interval in the argument of p exactly 

corresponds to the inequalities Q { °\ = Y 1 + \ti > 1 > Y 1 + = Q^'J- 

To estimate the supremum, we recall that p(dx) = (j) p (x)dx where <p p is bounded. Then 

K n {9,9')< ||^,||oo|0-0 , |&~ n+ \ nGN. (4.24) 

Now, by Corollary 3.8, Z n (0) < C$(p) n for some C < oo uniformly in p on convex sets N C 
{p € A 7 ^ : 3(p) < 3 c e} with finitely many extreme points and uniformly in 9 < 9q. Therefore, 
the right-hand side of (4.22) is bounded by B 00 {9q, \)\9 — 9'\ times a sum that converges when- 
ever 3 (p) < 3 c e, uniformly in p G A/", where AA is as above. This proves the desired claim. □ 

Let p 6 AT and let Qoo be the random variable defined in Section 3.3, independent of both 
the green sites and X a . Let us introduce the quantity 

B^(A)=E( J B 0O (Q oo ,A)), (4.25) 

The significance of B^ Q (X) is that it represents a stationary form of B 00 (-, A), i.e., B^ >0 (X) is a very 
good approximation of the probability P /9jA (^' ) n <3 = | a G ^ v) ), where 0' = Q { a X0+v) 
and where ^ 9 ^ is the quantity J^ 6 *') for trees rooted at a very far from 0. Let 

k p (\) =l([E( J B 00 (A + %,A))] 2 |g oo > l), (4.26) 

where X and Qoo are independent with distributions P = p and P, respectively. For critical 
distributions, B^ X} {\) and x p (A) are related as follows: 

Lemma 4.6 Le? p G A4 b be such that j(p) = 3 C . T/ien 

Proof. Since B^Q, A) — > as A j 0, we can expand on the right hand side of (4.19) to 
the second order of Taylor expansion, use that %{p) = P(Qoo > 1) and apply 63 (p) = 1 with 
the result 

BS (A) = A + (l-A)BSo(A)-^x p (A)[l + o(l)], A 1 0. (4.28) 
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(Here we noted that B OCl (X + ^Qoo) < B 00 {9b) allows us to estimate the error in the Taylor 
expansion by k p (X)B 00 (9i ) )0(1), which is x p (A)o(l) as A [ 0.) Subtracting (1 — \)B£ (\) from 
both sides and dividing by A, (4.27) follows. □ 

Note that, by the resulting expression (4.27), x p {\) / A tends to a definite limit as A J, 0. In the 
supercritical cases, on the other hand, Lemma 4.6 gets replaced by the following claim: 

Lemma 4.7 Let p,p' G A4 and define p a = (1 — a)p + ap'. Suppose that 3(p) = 3 C 
and 3(p a ) > 3c for all a G (0, 1]. Let -BJ-^O, a) denote the quantity -B^O) for the underly- 
ing measure p a . Then 

B » (0 -" ) = %3^ [1 + o<1)] - al °- (4 - 29) 

Proof. As in Lemma 4.6, we use that -Boo(#, 0, a) — > as a — > 0, where B^Q, 0, a) denotes the 
quantity B^Q, 0) for the underlying measure p Q . However, instead of (4.28), this time we get 

B^(0,a)(l-b}(p a ))=-^x Pa (0)[l + o(l)], a | 0, (4.30) 

where we again used that the error in the Taylor approximation can be bounded by x Pa o(l). 
Dividing by 3(p Q ) - 3c / 0, (4.29) follows. □ 



4.4 Critical asymptotics. 

The purpose of this section is to finally give the proof of Proposition 4.2. We begin by proving an 
appropriate upper bound on B^O, A). Note that, despite being used only marginally, equation 
(4.27) is a key ingredient of the proof. 

Lemma 4.8 Let p G M 9 satisfy $(p) = 3 C . For each 9 > 1 there is a K{9) G (0, oo) such that 

lim S up Boo( ^ A) <K{9). (4.31) 
A|o VA 

Proof. Let 3 = 3(p). We begin by proving (4.31) for 9 = 1. Let 

i{p) = %(p([l - -3-Qoo, 1]) 2 |Qoo > l) (4.32) 

and recall the definition of x p {\) in (4.26). Using the inequality B 00 (9, A) > B^l, A)1 {e>i} 
we derive x p (A) > i(p)B OQ (l, A) 2 . Inserting this in (4.27), we have 

B^o(X) < 1 — ~~^j~ L (p) ^oo(l, A) 2 [l + o(l)] , A 10. (4.33) 

Since the left-hand side is always non-negative, (4.31) for 9 = 1 follows with K(l)~ 2 = ^-/.(p). 

Next we will show that for any 9 < 9b, B OQ {9, A) is bounded above by a (^-dependent) multiple 
of -Boo(L A). Indeed, pick an e > such that 9^ — 9 > and let m be so large that (2.15) 
holds. Fix a directed path of m steps in T5 starting from the root. By conditioning on the event 
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that X a > x+ — e for all a ^ in the path, we have B^l, A) > p([x+ — e, l]) m B OQ (9, A), i.e., 

Boo(e, A) < C(0)Boo(l, A), 9 < 9 b , (4.34) 

with C{9) = p([x* - e, l])~ m < oo. 

As the third step we prove that (4.31) holds for values 9 in slight excess of 9 b . (The reader 
will notice slight similarities with the latter portion of the proof of Theorem 2.4(1).) Let e > 
be such that — e > 1 — By Corollary 2.5 and the fact that p £ we can assume 
that K e = p([x* - e, x*\) < 3. If 9 > 9 b is such that 9 e = x* - e + \9 < 9 b , then (4.19) and the 
bound &b(y) < by imply 

Boo(9, A) < A + (1 - A)6[k £j B 00 (^, A) + (1 - /c e )Boo(0 e , A)] , (4.35) 

because X + ^9 < 9 for all X in the support of p. Since (1 — X)bK e < bn e < 1, we have 

B.(^)< A+(1 - A ? (1 :;<y^ (1 - A) . <436) 

1 — (1 — X)bn e 

Dividing by ^A and taking A | 0, (4.31) follows with K{9) = 6(1 - K e )C(9 e )K(l)/(l - bn e ). 
Finally, it remains to prove (4.31) for general 9 > 9 b . But for that we just need to observe that 

Boo(9, A) < [1 - (1 - Xf +1 ] + (1 - X) bk+1 B 00 (9 b + 9b-\ X) (4.37) 

as follows by conditioning on the first k layers of T& to be green-free. By taking k large enough, 
9 b + 9b~ k is arbitrary close to 9 b , so the result follows by the preceding arguments. □ 

Lemma 4.8 allows us to write the following expression for l?oo(#, A): 

Lemma 4.9 Let p G M b satisfy l(p) = fc. Let e(A, 9) be defined by 

Boo(0, A) = ^j p (9)B^(X) + e(A, 9), (4.38) 

where ^> p (9) is as in (3.20). Then lim^o e(A, 9)X~ l l 2 = uniformly on compact sets of 9. 

Proof. Recall the notation Q n [ from (3.12), and let E„ denote the expectation with respect to 

(0) 

the measure P„ in (3.13). We will first show that 

B^e, A) = Z n (9)b n E^ (BooCQW, A)) + e n (A) (4.39) 

holds with an e n (A) satisfying lim^o £n{X)X~ l l 2 = for all n > 1. Let G n denote the ra- 
th generation of T&, i.e., G n = {a G T b : \a\ = n}, and let H n = U m <n ^m- Recall the 
notation for the object on the subtree of T b rooted at a and let be as described 
in (2.6). Given a a € G n , let ir(a) = {m k (a) : k = 1, . . . , n} be the path of connecting a to the 
root. 

A moment's thought reveals that, if f nH„ = (i.e., if there are no green sites in the first n — 1 
generations of T^), then in order for n / to occur, the following must hold: First, there 
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(0) 

is a a G G n , such that Q a , > 1 for all a' G vt((t). Second, the avalanche starting from this a 
with an initial amount reaches ^ . Introducing the event 

u n = [j n <s ± 0} n p| {qJ? > i}) , (4.40) 

we thus have 

W p ,x(Un) < #oo(M) < P ft A(W n )+P ftA ({^nE n + 0}). (4.41) 

Since P p (^nEI n / 0) = O(A), it clearly suffices to show that W\^(U n ) has the same asymptotics 
as claimed on the right-hand side of (4.39). 

Since U n is the union of b n events with the same probability, the upper bound 

P p ,a(Z4) < 6 n Z n (^)E(f)(S 00 (Q; t e) 1 , A)) (4.42) 
directly follows using the identity 

E p (£UQW,A) [J 1 {Q w> l} ) =Z n {9)K ( i\B 00 {Q%\)). (4.43) 

cr'eTr(o-) 

To derive the lower bound, we use the inclusion-exclusion formula. The exclusion term (i.e., the 
sum over intersections of pairs of events from the union in (4.40)) is estimated, using the bound 
in Lemma 4.8, to be less than K(9) 2 b 2n X, where = V 9 b . This proves (4.39). 

Since l{p)b = 1, Corollary 3.8 tells us that Z n (6)b n = ip p (9) + o(l). The final task is to 
show that (5^ (Qjf J, A)) can safely be replaced by its limiting version, B^X). We cannot 
use Corollary 3.6 directly, because 9 i— > i?oo(# 5 A) is known to be Lipschitz continuous only 
for 9 > 1. However, by Lemma 4.4 we know that i?oo(#, A) = A for 9 < 1, which means that we 
can write 

B 00 (9,X) = B 00 (9Vl,X)+[X-B 00 (l,Xp {e<1} . (4.44) 

Now, 5^,(0, A) = B^ V 1, A) is Lipschitz continuous in 9 for all 9 > 0, so by (4.21) and 
(3.18), 

EW(^(QS,A)) -i(^(Qoo,A))| < DB 00 (9,X)e-" n (4.45) 

where c > and D = D(9) < oo. To estimate the contribution of the second term in (4.44), 
we first note that A — i?oo(lj A) is a constant bounded between —B^f), A) and zero. Hence, we 
thus need to estimate the difference Fn\Qn \ < 1) — P(Qoo < !)• But that can be done using 
Proposition 3.5: Let k = [§J and use the monotonicity of 9 and (3.14) to estimate 

K 9) (Q!3 < 1) - HQoo < 1)| < {Qt\ > 1) - P(Qg! > 1) < A"e-« n ~ k \ (4.46) 

where A" = A/(l — e"^). By combining all the previous estimates and invoking (4.31), we find 
that the difference (B^iQ^, A)) - B^(X) is proportional to e" f ' n \/A, where q' > 0. Using 
this back in (4.39) the claim follows by taking the limits A j and n — > oo. □ 

Lemmas 4.8 and 4.9 finally allow us to prove Proposition 4.2: 
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Proof of Proposition 4.2. Note that, by using (4.38) in (4.26) and the definition of c p in (4.11), 
we have 

b -Z±„ p (\) = B* 00 (\) 2 cf + o(\), A |0. (4.47) 

Then the fact that B^(X) tends to zero as A j forces, in light of (4.27), that ^x p (A) -> 1 
as A J, 0. This in turn gives that 

B£,(A) = VA(c p + o(l)), AjO. (4.48) 
Plugging this back in (4.38) proves the desired claim. □ 

4.5 Supercritical case. 

Here we will indicate the changes to the arguments from the previous two sections that are needed 
to prove Proposition 4.3. We begin with an analogue of Lemma 4.8: 

Lemma 4.10 Let p,p' £ and define p a = (1 — a)p + ap'. Suppose that $(p) = $ c 
and }{p a ) > Icfor all a 6 (0, 1]. Then for each 9 > 1, there is a constant K'(9) 6 (0, oo) 
such that 

limsup^f'"' 1 '^^^). (4.49) 

alO KPa) - 3c 

Proof. The only important change compared to the proof of Lemma 4.8 is the derivation of 
the bound for 6=1. Indeed, in this case we use that x Pa (0) > B^O, a)-Boo(l, 0, a) in (4.29), 
where B^l, 0, a) is the quantity B^O, A) for A = 0, 9 = 1 and p = p a . Applying -B^(0, a) > 
for all a G (0, 1], as follows by Theorem 3.1(2), we find that (4.49) holds with if' (1) = 
Once we set A = 0, the rest of the proof can literally be copied. □ 

Next we need to state the appropriate version of Lemma 4.9: 

Lemma 4.11 Let p,p' G Ai^ and define p a = (1 — a)p + ap'. Suppose that %{p) = 3 C 
and l{p a ) > 3c far alia G (0, 1]. Then 

yg^ .^ Mj^l ^,), ai0 , (4.50) 

3(Pa)-3c 3(Pa)-3c 
where K a is the expectation corresponding to Ffor measure p a . 

Proof. Also in this case the required changes are only minuscule. First, we have an analogue 
of (4.39), 

¥ Pa (\^\ = oo) = b n Z^\9)E^ a (P, a (|^0| = oo))+e' n (a), (4.51) 

where En} a is the expectation and Zn"* the object Z n {9) for the underlying measure p a and 
where e' n {a) is the quantity in (4.39) for A = and p = p a . We claim that 

lim e ' n [ a ^ = (4.52) 
<4o i(p a ) - 3 C 
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for all finite n > 1. Indeed, the entire derivation (4.40-4.45) carries over, provided we set A = 
0. The role of the "small parameter" is now taken over by l(p a ) — 3c- A computation shows 
that e„(a) = 0((3(Pa) - 3 C ) 2 ) as a j 0, proving (4.52). 

To finish the proof, it now remains to note that b n Z^ a] {9) -> b n Z^\o) as a [ and that, by 
Corollary 3.8 and the fact that j(p) = 3 C , we have b n zi p) (9) = ip p (9) + o(l) as n — > oo. □ 

Recall the definition of c p in (4.11). To prove Proposition 4.3, we will need to know some 
basic continuity properties of c p in p. Note that these do not follow simply from the continuity 
of oi i ^ Vp a (0)> because also the expectation E in (4.1 1) depends on the underlying measure. 

Lemma 4.12 Let p,p' £ M!" be such that p a = (1 — a)p + ap' satisfies l{p a ) > for 
all a £ [0, 1]. Let c p be as in (4.11). Then lim^o c Pa = c p . 

Proof. Let ip* a (9) = K Pa (ip pa (X + \0)). In general, ip Pa (6) is Lipschitz continuous for 9 > 1. 
Thus, ip Pa converges uniformly to tf) p on compact sets of 9. Hence, we just need to show 

HmE a (^(Qoo) 2 |Qoo > 1) = E(^(Qoo) 2 |Qoo > 1). (4.53) 

Choose n > 1 and replace K a , E and Qoo by their finite-n versions. By Corollary 3.6, the error 
thus incurred is uniformly small in a G [0, 1]. Hence, it is enough to show that 

W(oa) 2 |oiS > i) = W(qJ?i) 2 |o21 > i), (4-54) 

for some G [1,^&], where E„ ; p denotes the expectation with respect to P„ for measure p. 
However, in (4.54) only a finite number of coordinates are involved and the result follows. □ 

With Lemmas 4. 10, 4. 12 and 4. 1 1, we can finish the proof of Proposition 4.3: 
Proof of Proposition 4.3. From (4.50) we have 

I '" 1 x Pct (0) = ^(0,a) 2 C ; ( 2 + O (3(p a )-3 c ), a 10. (4.55) 

12, we have 

6c 2 + o(l), a|0. (4.56) 



2 

Using this in (4.29) and invoking Lemma 4.12, we have 



3 (Pa) ~3c 

The proof is finished by plugging this back into (4.50) and invoking the continuity of a ^ 



5. Coupling argument 

5.1 Coupling measure. 

The goal of this section is to define a coupling of the measures and P^f ^ that appear in (3. 14). 
As the first step, we will write Fn \-) as the distribution of a time-inhomogeneous process. To 
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have the process running in forward time direction, we will need to express all quantities in terms 
of the (more or less) original variables (Xk), which relate to the F's through 

X k = Y n _ k+1 or Y k = X n _ k+1 , 1 < k < n, (5.1) 

(0) 

see Section 3.3. Abusing the notation slightly, P„ '(■) will temporarily be used to denote the 
distribution of the Xi, . . . , X n as well. We will return to the V s in the proofs of Propositions 3.5 
and 3.7. 

Let Z n (9) be as in (2.8) and note that, since p G , we have Z n {9) > for all n > and 
all 9 > 1. Given 1 < k < n — 1 and, for k > 1, a sequence (Xl, . . . ,X k _{) G [0, 

let C(-) = O ■ l^i. • • • . be § iven b y 

e 



where the indicator ensures that we are not dividing by zero. The (X\, . . . , X k _{) -dependence 
of t^ k will be often left implicit. 

To interpret these objects, let us consider the case k = 1. Suppose that we wish to elucidate the 
distribution of X\ knowing that the process will survive long enough to produce an X n -\. (The 
variable X n corresponds to Y\, which will be uncorrected with the other Y's.) The only prior 
history we know is the value of 9; obviously we are only interested in the case 9 > 1. The total 
weight of all configurations is just Z n _\{9); hence the denominator of (5.2). Now, if X\ takes 
value x, the weight of configurations in which the process survives is like the weight of a string of 
length n — 2 with an effective "9" given by x + ^9. Hence Z n -2 {x + \9) in the numerator. (Notice 

thatifx+^ < 1, this automatically vanishes.) We conclude that P^ {X\ G Ax) = tf\{x)p(Ax). 

A similar reasoning shows that the probability of {X k 6 dx} given the values of Xi, . . . , X k ~i 
equals t^\(x)p(dx). This allows us to view P^ as the distribution of an inhomogeneous process: 

Lemma 5.1 For all 9 > 1, all n > 1 and all Borel-measurable sets 4 C [0, 1]", 

n— 1 

PW(A) = E(1 a II €!k^k\Xu ■ ■ .,**_!)) • (5-3) 
fc=i 

Proof. The result immediately follows from the formula 

n— 1 ^ , n— 1 x 

now.-.*-.) = n w<« l2 i> )• < 5 - 4 > 

the identity Qjjf ) = X k + |qJJi and the definition of P^ e) (•), see (3. 13). □ 

Next we will define the coupled measure. The idea is to use the so-called Vasershtein coupling, 
see [8], which generates new (coupled) pairs from the "maximal overlap" of the individual dis- 
tributions. Let 9, 9' > 1, and suppose that the corresponding sequences X = (X\, . . . , X k _i) € 
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[0, and X' = (X[, . . . ,X' k _ x ) £ [0, have been generated. Assume also that a se- 
quence (loi, . . . ,Uk-i) € {0, satisfying uji < ~\{x e =x' e } for all 1 < ^ < fc — 1 has been 
generated. (This sequence marks down when X^ was coupled with X' e . Note that we could 
have that X? = X' e even when Xp and X' £ are not coupled.) Let t be the quantity for the 
sequence X and let t' be the corresponding quantity for the sequence X' . Let 

R( ■ ) = R^t\' l^i> • • • i x k-i', X[, ■ ■ ■ ,X' k _i,uJi, . . . ,cok-i) (5-5) 

be the transition kernel of the joint process, which is a probability measure on [0, 1] x [0, 1] x {0, 1} 
defined by the expression 

/ / -,\ j t(x) At'(x) p(dx)8 x (dx'), ifu = l, 

R(dxxdx x M ) = j^ [t(x) _ nx)]+[n ^_^ )]+p(da;)p(dx0) ifw = _ (5-6) 

Here t(x) A t'(x) denotes the minimum of t(x) and t'(x) and [t(x) — t'(x)] + denotes the positive 
part of t(x) — t'(x). The quantity q = q^f-x X' * s gi ven by 

, = / tW A^)^) = l-/[ tW -^)] + ^). (5.7) 

The interpretation of (5.6) is simple: In order to sample a new triple (X^, X' k ,Uk), we first 
choose Wjfc G {0, 1} with Prob(a;fc = 1) = q. If w& = 1, the pair (Xk,X' k ) is sampled from 
distribution ^t(x) A t'(x) p(dx)S x (dx') — and, in particular, gets glued together with X' k — 
while for the case u>k = we use the distribution in the second line of (5.6). 

Remark 5. It turns out that whenever the above processes X and X' have glued together, they 
have a tendency to stay glued. However, the above coupling is not monotone, because the pro- 
cesses may come apart no matter how long they have been glued together. Our strategy lies in 
showing that q tends to one rapidly enough so that the number of "unglueing" instances is finite 
almost surely. 

(9 6') 

Let Pn ' (■) be the probability measure on [0, 1]™ x [0, 1]™ x {0, 1}™ assigning mass 

„ rt— 1 

P<W(B) = J2 P (dx n )p(dx' r ^ {u}n=1} J] <f ; 2,^dz fc xd4x{u; fe }) (5.8) 
(«*) *=i 

to any Borel-measurable set B C [0, 1]™ x [0, 1]™ x {0, 1}™. Here R^Q^ (dx h xte k x.{wk}) = 

(0 6 r ) 

R n ' k (dxk x dx f k x {u;/J|xi, . . . , x^-i] x[, • • • , x^^Ui, . . . , Wk-i)- As can be expected from 

(0) (6 r ) (6 6 f ) 

the construction, P n (*) and P n } (•) are the first and second marginals of P n ' 1 (•), respectively: 
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Lemma 5.2 Let 9, 9' > 1. Then 

VP' 9 "* (A x [0, 1]" x {0, 1}") = PW (A) (5.9) 

and 

Pif' e ">( [0, 1]™ x A x {0, 1}") = Pf (5.10) 
for all Borel-measurable A C [0, l] n . 

Proo/ To prove formula (5.9), let X = {X 1 ,..., X k -i) and X' = (X[ , . . . , X'^) be two 
sequences from [0, l]* -1 . If > 1 and the same holds for the corresponding quantity for 

the sequence X', let t(-) = t^ k (-), t'(-) = t^2('). and let R ( m ) and 1 be as in ( 5 - 6 ) and ( 5 - 7 )' 
respectively. Using (5.7) we have, for all Borel sets C C [0, 1], 

V / R(dxxdx' x{oj}) = / (t(x) At'(x) + [t(x) - t' (x)} + ) p(dx) = / t(x)p(dx). 

(5.11) 

In other words, the first marginal of the coupled process is a process on [0, 1] with the transition 

(0) 

kernel t(-)p(-), which, as shown in Lemma 5.1, generates P„ ; . This proves (5.9); the proof of 
(5.10) is analogous. □ 

Clearly, the number q represents the probability that the two processes get coupled. The fol- 
lowing lemma provides a bound that will be useful in controlling q: 

Lemma 5.3 Let 9, 9' > 1, 1 < k < n - 1 and X = (X u . . . , X k -i) € [0, l]^ 1 and X' = 
(X[, . . . , X' k _]) G [0, Let Q be the quantity corresponding to X and let Q' be the 

quantity Q k _i corresponding to X'. If Q A Q' > 1, then 

(6,V) Z n _ k (QAQ>) 
q n,k-,x,x>> Zn _ k{QvQ/ y ^- 12 -> 

(0) 

Proof. Let t be the quantity t n k for the sequence X and let t' be the corresponding quantity for 
the sequence X'. By inspection of (5.2) and monotonicity of 9 Z n {9), 

and similarly for t'(x). From here the claim follows by integrating with respect to p(dx). □ 
5.2 Domination by a discrete process. 

The goal of this section is to show that the coupled measure defined in the previous section has the 
desirable property that, after a finite number of steps, the processes X and X' get stuck forever. 
Since the information about coalescence of X and X' is encoded into the sequence uj, we just 
need to show that, eventually, u) k = 1. For technical reasons, we will concentrate from the start 
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on infinite sequences (oJ k ) ke ?$: Let Pn' 6 \-) be the law of (uJk)keN £ {0, 1} N induced by the 
distribution Fn' e \-) and the requirement Pn' 9 \to k = 1, k > n) = 1. 

The coalescence of X and X' will be shown by a comparison with a simpler stochastic process 
on {0, 1} whose law will be distributionally lower than Pn ' (•), i.e., in the FKG sense. Let =<! 
be the partial order on uj,uj' G {0, 1} n defined by 

co 4 J 44> uj k <uj' k , k>l. (5.14) 

Next, note that, by x+ > ^p, we have 1 — 6(1 — x*) > 9b — 1. Choose a number 5 P G 
(Ob — 1, 1— 6(1— x+)) and, noting that p([l— ^-jp-, £*]) > 0, define a collection of weights (A p (s)) 

V s ) fc > s d-P^b-" p([l ~ V' 1 *]) 

Note that s A(s) is increasing. It is also easy to verify that A p (-) G (0, 1], so any of these 
weights can be interpreted as a probability. This allows us to define a process on (oJ k ) k€ -^ €= 
{0, 1} N , with the transition kernel 

p p (u' k = 1 \lo[, . . .,<J k _ x ) = A p (min{0 < j < k - 1: J k _^_ x = 0}), (5.16) 

where, for definiteness, we set uo' = 0. Let P p (-) denote the law of the entire process with 
transition probabilities p p ( ■ \ ■ ) and "initial" value oj' = 0. 

Proposition 5.4 Let p G and let 5 P be as above. For all n > 1 and all 9, 9' with 1 < 9, 9' < 

(9 9') ~ 
9b, the measure P n ' (•) stochastically dominates P p (-) in partial order ^. 

Let S p be fixed for the rest of this Subsection. In order to give a proof of Proposition 5.4, we 
first establish a few simple bounds. 

Lemma 5.5 Let p G Ai^ and let 5 P be as above. Let n > and suppose 9,9' > 1 satisfy < 



-9' < 5 p b~ k for some k > 0. Then 



Proof. Consider a configuration X±, . . . , X n which contributes to Z n (9) but not to Z n (9'). This 

implies that there is an £ G {1, . . . , n} where > 1 but Q\ < 1. With this in mind, we 
claim the identity 

n n n r i i r ra 

n \q^>i } - n i { q^)>i}=e n i w ^)>i} vr^^} n \ Q ^>i } 

m=l m=l t=\ m=l J L m=£+l 

(5.18) 

Thence, 

n , e-i \ 

Z B (0) - Z B (^) = 1 wf )<1 < Q (.) } IT 1 ^')>i } j- ^ 19 ) 

1=1 ^ m=l ' 
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Since 6-6' < 5 p b~ k , we have qP - 1 < qP - Qf' ] < 5 p b~ k ~ e for any i contributing on the 

right-hand side. In particular, we have qP < 1 + ^, which implies Z n -i(QP) < Z n _£(l + ^). 
Then 

n , l-\ 

z^-z^^Zn^ + ^Hp^-iQ^-lQtl)) II 1 

i=l ^ m=l 

or, replacing p([l — ^Q\_ x , 1 — ^Q\_[)) by its maximal value, 



«sr>>i} )■ (5 - 20) 



Z n {6) - Z n {6>) < V Z n _ £ (l + ^(0') sup p([l - f , 1 - £ )). (5.21) 

On the other hand, by simply demanding that X^> 1 -g-^ (which implies Q e > 1 + in 

(2.10) we have for all 1 < £ < n that 

Z n (0') > Z n _ e (l + £)p([l - ^XiDZe-iiO'). (5.22) 
Using (5.22) in (5.21), and applying (5.15), we have 

Z n {6) - Z n {6>) < l ~ X j k) Z n (6% (5.23) 

whereby the claim directly follows. □ 
Next we prove a bound between kernels (5.6) and (5.16): 

Lemma 5.6 Let 1 < k < n-\andletJ = (uj[, . . . ,to' k _ 1 ) G {0, X = (X 1 , . . . ,X fc _i) G 

[0,l] k -\ X' = (X[,...,X' k _ 1 ) G [0, l]*" 1 and w = (wi, . . . , w fc _i) G {0, l}*" 1 . For 

all 6,6' > 1 araci a// £ = 1, . . . , k — 1, let Q^' correspond to X via (2.5), a«<i Ze? cor- 
respond to X'. Suppose that 

Q ( p > 1, Qf } > 1 W ^ < wj < 1 {Xj . =x j } , i = 1, ...,/«- 1. (5.24) 
fjf -R^ j^/ w quantity defined in (5.8), ?/z£?2 

= 1}) > JvK = 1 K, • • • (5-25) 

for all 6,6' with l<6,6' <6 b . 

Proof. Note that, since 1 < 6,6' < 6 b and 1 + S p > 6 b , we have 1 < Qf ] ,Q { p < 1 + 5 P and 
thus — Q\ | < (5 P for all i = 1, . . . , k — 1. This allows us to define the quantity 

s = max{f: < £ < k, \Q^ - Qf\\ < Spb'^}. (5.26) 

By Lemmas 5.3 and 5.5, we have R({uJk = 1}) > \>( s )> where R(-) stands for the quantity on 
the left-hand side of (5.25). Recall our convention uj' = and let 

s' = min{0 < j < k - 1 : <4-j-i = 0}. (5.27) 
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In other words, s' is the length of the largest contingent block of l's in u' directly preceding u>' k . 

We claim that s > s'. Indeed, by our previous reasoning, \Q k _ s i_i — Qj.-s'-il — *V ^ our 
assumptions, 1 = u'j < 1{x,=x'.} an d, therefore, Xj = X' 3 for all j = k — s', . . . , k — 1. 
This implies 

iQiV^Pil (5-28) 
and hence s > s'. Using that s' is the argument of A in (5.16) we have R{{uOk = 1}) > -V( s ) > 
Xp(s') = p p ( uj' k = l\u>[, ... , u k _i)- This proves the claim. □ 

Now we are ready to prove Proposition 5.4: 

Proof of Proposition 5.4. The inequality (5.25) is a sufficient condition for the existence of so- 
called Strassen's coupling, see [8]. In particular, the inhomogeneous-time process generating the 
triples (X k ,X' k ,Lo k ) can be coupled with the process generating u> k in such a way that (5.24) holds 
at all times less than n. The (to, uj') marginal of this process will be, by definition, concentrated 

(6 6') 

on {oj ^ lo'}. Since oo k = 1 for k > n, Pn ' -almost surely, the required stochastic domination 
follows. □ 

5.3 Existence of the limiting measure. 

The goal of this section is to show that, under proper conditions, the process oo' with distribu- 
tion Pp(-) equals one except at a finite number of sites. Then we will give the proof of Proposi- 
tion 3.5. Let 

_/(l-A p (n))nj^A p (fc), ifnGNU{0}, 
Ulfc=o A p(*0> ifn = oo, 

and observe that p n is the probability of seeing a block of l's of length n in the prime configura- 
tion. We begin with an estimate of \{k): 

Lemma 5.7 For each p € M!", there is C{p) < oo and w > such that 

1 - \ p {k) < C(p)e~™ k . (5.30) 

Moreover, the quantity C(p) is bounded away from infinity uniformly in any subset M C 
with finitely many extreme points. 

Proof. Let (f) p be the density of p with respect to the Lebesgue measure on [0, 1]. Then 

sup p{[l - f , 1 - % )) < <yrl0joc. (5.31) 

e-e'<s p b- n 

The claim then follows by inspection of (5.15) with w = log b and an appropriate choice of C (p). 
The bound on C(p) is uniform in any N with the above properties, because the bound \\4>p\\ p < oo 
is itself uniform. □ 

The preceding estimate demonstrates that the discrete process locks, and in fact does so fairly 
rapidly. Indeed, we now have p^ > 0, which ensures that eventually the configuration is all ones, 
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and further that the p n tend to zero exponentially. It remains to show that the waiting times till 
locking are themselves exponential. 

Lemma 5.8 Let p G M h and, for n > 1, let £(n) = {J G {0, 1} N : Uj = 1, j > n}. 
Let ao > be such that <p(a) = ^ <fe<oo e a ^ k+1 ^Pk < oo for all a G (0, ao). Then 

P p (£(n) c ) <ne^, n > 1, (5.32) 

//(p) = S up{a > 0: <p{a) < l}. (5.33) 

We note that both quantities ao and p,(p) wee nontrivial. Indeed, ao > w > and, since 
can be written as p^ = 1 — J2 n >oP n > 0> we nave tnat M/ ) > 0- 

Proof. An inspection of (5.16) shows that "blocks of l's" form a renewal process. Indeed, sup- 
pose £g for £ = 1, . . . , k — 1 mark down the lengths of first k — 1 "blocks of l's" including the 
terminating zero (i.e., & = n refers to a block of n — 1 ones and followed by a zero). Denot- 
ing Nk-i = YljZi Cj> tne ^-th block's length is then 

C fc = min{j > : ^ +Wfc i = 0}. (5.34) 

As is seen from (5.16), (&) can be continued into an infinite sequence of i.i.d. random variables 
on N U {oo} with distribution Prob(^ = n + 1) = p n , where p n is as in (5.29). The physical 
sequence terminates after the first = oo is encountered. Let Q n {k) be the event that £ij • • • > £fc 
are all finite and Yli=i & > n - Then, clearly, £(n) c = Ufc=i S n (k). 

The probability of £/„(£;) is easily bounded using the exponential Chebyshev inequality: 

Prob(g n (fc)) < cp(a) k e- an , 0<a<a . (5.35) 

Noting that Ylk=i ( / 7 ( Q; ) fc < ^ for a < p(p), the claim follows. □ 
Now we are finally ready to prove Proposition 3.5: 

Proof of Proposition 3.5. Let p G and n be fixed. Let k < n and suppose that / is a function 
that depends only on the first k of the Y -coordinates. Let 9o > @b and let 9, 9' G [1, 9q]. Noting 
that p£\|Q$n G dQ) = P { ® m (-), we have 

KUf)-^'\f)\ < <l 1 (\^n (:k "\f)-H?'Hf)\). (5.36) 

Since Q n \_i n G [1, #o] by our choice of 0, we just need to estimate |E„ (/) — E„ (/) | by the 
right-hand side of (3.14) for all 9, 9' G [1, 9 ]. 
Introduce the quantity 

D n (/) = sup{ |E(f ) (/) - Ef ) (/) | : 9, 9' G [1, 9 ] } . (5.37) 

We need to show D n (f) is exponentially small in n. By Lemmas 5.1, 5.2, and Proposition 5.4, 
the probability that X, / X- for some n — k < i < n under the coupling measure P„ ' (•) is 
dominated by the probability that uj^ = for some n — k < i < n under P p (-). Since / depends 
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only on the first k of the Y variables (i.e., the last k of the X variables), the coupling inequality 
gives us 

|E(?) (/)-Ef)(/) | <2\\f\\ 00 P p (£{n-kf), (5.38) 

where £(n — k) is as in Lemma 5.8. 

Let \i = p,(p) be as in Lemma 5.8. Then (5.32) and (5.38) give 

D n (f) < 2||/||oo(n - k) e"^- fc ) < 4( M e)- 1 ||/|| 0O e-^(n-*) j ( 5 .39) 

This proves (3.14) with ( = ^p and A = 4( / ue)~ 1 . The bounds ( > and A < oo are 
uniform in sets M C with finitely-many extreme points, because the bound p(p) > is itself 
uniform. The existence of the limit (3.15) and its independence of is then a direct consequence 
of (3.14). ' □ 

5.4 Distributional identity. 

Here we will show the validity of the distributional identity (3.19). The proof we follow requires 
establishing that the distribution of Qoo has no atom at = 1: 

Lemma 5.9 Let p G M b . Then P(Qoo = 1) = 0. 

Proof. Notice that the almost-sure bound < Qoo < Q^i holds for all n > 1, with T 

Qoo and Q<M [ Q oo as n -> oo. Therefore, 

HQoo = 1) = lim P(QW < 1, > 1). (5.40) 

But Y\ is unconstrained under P(-) which by < Q^ b \ — Qn \ ^ (@b — l)b~ n allows us to write 

P(QS < 1, Q ( n"i > 1) < 1-h-s. of (5.31). (5.41) 
Hence, ¥(Q^ < 1, > 1) -> as n -> oo and we have P(Qoo = 1) = 0, as claimed. □ 

Proof of Proposition 3.7. Let X be a random variable with distribution P(-) = p(-), independent 
of Yi, Y2, ■ ■ ■ , and let 6 > 1. For all a G P, define the (distribution) functions 

Fi e) (a)=^(Q ( V>a). (5.42) 

and 

(a) = P ® pW (x + ^ > a, Q $ > 1 J . (5.43) 

Since \ = Q^+i 2> = ^i and ^i + i^i+i 2 = Qrf+i i> tnese functions obey the relation 

FW(a) = Fi 9 \l) F^a), n > 1, a G R. (5.44) 
Let F(a) = P(Qoo > a) and let 

F(a) = P®P(X + %^ > a, Qoo > 1 ). (5.45) 
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Both F(-) and F(-) are non-increasing, left-continuous and they both have a right-limit at ev- 
ery o £ 1. In particular, both functions are determined by their restriction to any dense subset 
of R. The proof then boils down to showing that there is a set A C R dense in R such that 

lim F^\a) = F(a) a G A U {1}, (5.46) 

and 

lim Ff)(o) = F(a), a G A. (5.47) 

Indeed, then (5.44) implies F(a) = F(l)F(a) for all a G A which by continuity extends to 
alia G R, proving (3.19). 

Let A be the set of continuity points of both F(-) and F(-). Clearly, A° is countable and 
hence A is dense in R. The limits in (5.46) will be taken in too stages; first we take the limit of 
the distribution and then that of the event. Since < \ < Q^i f° r an y rn < n,we have, 
by (3.15), 

HQm,i >a)< liminf F^{a) < lim sup F^ (a) < F(Q% b } > a ) ( 5 -48) 

n >oo n — >co 

for all 6 > 1 and all m > 1. The m —> oo of the extremes exists by monotonicity. Since 
Qmi > <3oo» the right-hand side converges to F(a). As for the left-hand side, it is clear that the 
event {Qoo > a} implies that, eventually, > a} occurs. Thus the limit of the extreme left 

is at least as big as F(Q OQ > a). However, the latter equals F(a) because, by assumption, a is 
a continuity point of F. This proves (5.46). The argument for the limit (5.47) is fairly similar; 
the right-hand side will directly converge to F(a), while the limit of the left hand side will be no 
smaller than P <g> F(X + -^Qoo > a, Qoo > 1). However, by Lemma 5.9 we have that P(Qoo = 
1) = and thus the limit equals F(a), because a G A. □ 

Proof of Corollary 3.6. The proof of Q\ [ — ► Qoo is immediate from (5.46). To prove (3.18), 
we note that (3.17) and (2.5) imply the deterministic bounds 

< C||/|U6-^o, (5-49) 

and 

|/(Qoo) - f{Q { ni)\ < C||/||oo6- n (9o, (5-50) 
where we used that x < Oo for 6 < 9q. The bound (5.49) implies that 

- <(f(Q^))\ < C'WfWooe-^, (5.51) 

where C' < oo and rj > 0, while the bound (5.50) guarantees that E(f(Q ao )) can be replaced 
by E(f(Q^\')) with a similar error. Then (3.18) with 2n replacing n boils down to the estimate 
of 

ES(/(QS))-E(/(QS)) . (5.52) 
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But, by Proposition 3.5, the latter is bounded by A||/|| 00 e _ '> ri . Combining all of the previous 
estimates, the claim follows. □ 

Proof of Corollary 3.8. We begin by showing that l(p) = P(Qoo > 1)- Indeed, we can use 
that Z n (9) = for 9 < 1 to compute 

E(Z n (Qoo)) =P(Qoo > l)E®E(z n -!(X + fan) \q oo > l) 

= P(Qoo > 1) I(Z n _ 1 (Q oc )) = • • • = PtQoo > l) n+1 , 

where we used Proposition 3.7 to derive the second equality. From here }{p) = P(Qoo > 1) 
follows by noting that P(Qoo > 1)Z„(1) < E(Z n (Q 00 )) < Z n (6 b ) and applying Theorem 2.4(1). 
In order to prove the existence of the limit (3.20), we first notice that 

^ =*&(<&,.>!)■ <5-54) 

Next we claim that P„+i(Q„+i i > 1) — l{p), for 9 > 1, decays exponentially with n. Indeed, 

let # > Ob and G [1, 9 ], pick fe = [§J , use Q^J < Qjfj^ < and apply Proposition 3.5, 
to get 

P(Qg > 1) - Ae~< k < ^UQ^ > 1) < P(Qg } > 1) + Ar<* (5.55) 

where ^4 < oo is proportional to ^4(p, 0q) fr° m (3.14). On the other hand, we clearly have 

P(QS > l) < P(Qoo > 1) < HQ k y > 1) • (5.56) 

But the right and left-hand sides of this inequality differ only by F(Q^\ < 1, > 1), which 

can be estimated as in (5.41) by a number tending to zero exponentially fast as k — > oo. From 
here we have 

where A' = A'(p, 9q) < oo and (' = ('(p) > 0. The uniformity of these estimates is a conse- 
quence of the uniformity of the bounds A < oo and ( > and that as in (5.41). 

The existence of the limit (3.20) for 9 G [1, 9q] is a direct consequence of (5.57) and the identity 

M 9) = M Zn (Wr = }nn \[ = fi (5.58) 



<A' e -f ,n , 9e[l,fio], (5-57) 



and the fact that the corresponding infinite product converges. For 9 < 1 we have Z n {9) = and 
the limit exists trivially. To prove that 9 ^ ^ P {9) is Lipschitz continuous for 9 > 1, we first note 
that, by (5.23) and the result of Lemma 5.7, 

\Z n {9) - Z n {9')\ <C\9-9'\iP p (9v)l{ P y n , 9,9'e[l,9 Q ], (5.59) 

where C = C(p, 9q) < oo is on sets M C with finitely many extreme points. From here the 
bound in part (2) directly follows. 
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Let Z„ (9) denote explicitly that Z n {9) is computed using the underlying measure p. The 
continuity of a ^ ip Pa {9) then follows using three facts: First, a i-> Z^ pa \9), being an expecta- 
tion with respect to p™, is continuous. Second, by Theorem 2.4(2), a i— > is also continuous. 
Third, the infinite product (5.58) converges uniformly in a. □ 
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